Two-loop approximation in the Coulomb blockade problem 
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We study Coulomb blockade (CB) oscillations in the thermodynamics of a metallic grain which 
is connected to a lead by a tunneling contact with a large conductance go in a wide temperature 
range, Ecgoe~ 90 ^ 2 < T < Ec, where Ec is the charging energy. Using the instanton analysis and 
the renormalization group we obtain the temperature dependence of the amplitude of CB oscillations 
which differs from the previously obtained results. Assuming that at T < Ecgte" 90 ^ 2 the oscillation 
amplitude weakly depends on temperature we estimate the magnitude of CB oscillations in the 
ground state energy as Ecgt e ~ 90 ^ 2 ■ 



PACS numbers: 73.23Hk, 73.50Bk, 73. 21. La 

I. INTRODUCTION AND MAIN RESULTS 

The study of electron-electron interactions in meso- 
scopic systems has been at the focus of experimental 
and theoretical interest over the past two decades. One 
of the most striking consequences of electron interac- 
tions at low temperatures is the phenomenon of Coulomb 
blockaded For example, thermodynamic quantities of a 
metallic grain which is connected by a tunneling contact 
to a metallic lead and is capacitively coupled to a gate 
exhibit oscillatory dependence on the gate voltage. This 
can be observed experimentally, by measuring the differ- 
ential capacitance of the grain.el 

At low temperatures the Coulomb interaction of elec- 
trons in the grain can be described within the framework 
of the constant interaction model 



He — Er 
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(1) 



where Ec is the capacitive charging energy of the grain, 
N is the operator of the number of electrons in it, and q 
is the dimensionless gate voltage. 

For a grain with a vanishing mean level spacing and 
for a tunneling contact with a large dimensionless con- 
ductance go — 2nh/e 2 R 3> 1, where R is the resistance of 
the contact, the amplitude of Coulomb blockade oscilla=j 
tions in the ground state energy is exponentially small,™ 
cx Ec exp(— The magnitude and temperature depen- 
dence of Coulomb blockade oscillations in this problem 
were studied in many theoretical Wftrks using the renor^ 
malizatioflr-group (RG) treatmentDB and the instantontl 
approachQH to the dissipative action of Ref. ^. 

In the instanton approach the cost of one instanton is 



tx exp(-^p), where g(T) « g -*^ 
malized conductance. The instanton gas can be consid- 
ered as non- interacting at relatively high temperatures, 
when the renormalized conductance g(T) is still large. In 
this regime the thermodynamic potential depends sinu- 



-2 In 



gpE c U 



is the renor- 



soidally on the gate voltage 



m 



n osc (q) = -E C {T) cos(27rg), 



(2) 



where Ec (T) is the renormalized temperature-dependent 
charging energy. 

Integration over the massive fluctuations around the 
instanton in the Gaussian approximation leads to the log- 
arithmic temperature dependences of the renormalized 
charging energy at relatively high temperature. This re- 
sults in the estimate i?c5o ex P(~^f) f° r the amplitude of 
the Coulomb blockade oscillations in ground state energy. 

In this paper we evaluate non-Gaussian corrections for 
the fluctuations about the instanton configurations in the 
lowest order of perturbation theory in l/go- We show 
that these corrections diverge at T — > and significantly 
modify the magnitude and the temperature dependence 
of the preexponential factor of the Coulomb blockade 
oscillations even at relatively high temperatures, where 
the renormalized conductance g(T) is large and non- 
interacting instanton gas approximation is still valid. We 
then apply a two-loop renormalization group (RG) to 
determine the preexponential factor in the renormalized 
charging energy beyond the region of applicability of per- 
turbation theory. 

The main result of this paper is the following expres- 
sion for the renormalized charging energy in Eq. (0), 
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Here g(T) is the renormalized temperature dependent 
conductance, and Xn = 4tt< With sufficient accu- 
racy g(To) can be found using the perturbation theory, 
g(To) =go — 2 In go- At lower temperatures, but still such 
that g(T) 3> 1 the renormalized conductance g(T) can be 
found from the implicit relation, 



g(T)= g Q -2 In 
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(3b) 
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which solves the two-loop renormalization group equa- 
tion.! 

Equations (||) and (||) rely on the non-interacting 
instanton gas approximation and are valid at T > 
Ec9o e ~ 90 ^ 2 ■ Assuming that at lower temperatures 
the amplitude of Coulomb blockade oscillations is only 
weakly temperature-dependent we obtain the estimate 
for the magnitude of oscillations in the ground state en- 
ergy, Ecgoe~ 9 "/ 2 . This is by a factor of go greater than 
the result of Ref. ||. 

The paper is organized as follows: In Sec. [n] we de- 
scribe the dissipative action approach to the problem. In 
Sec. p| w e discuss the instanton gas approximation. In 
Sec. f II Aj we treat the fluctuations around the instan- 
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tons in the Gaussian approximation, and in Sec 
we obtain the leading 1/go correction to the Gaussian 
result. In Sec. |y| we use the two-loop renorm alizati on 



group to extend the perturbative results of Sec. pi B 
the low temperature regime hi(Ec/T) ~ go 
are summarized in Sec. [v| 



to 



Our results 



II. 



PARTITION FUNCTION 



Following Ref. g the partition function of the system 
can be written as an imaginary time functional integral 
over the auxiliary field 4>{t) which decouples the Coulomb 
interaction Eq. (Jjj). Different configurations of the field 
4>(t) fall into different topological classes labeled by the 
winding number W, characterizing the boundary condi- 
tions in imaginary time: 4>w( T + 0) = 4>w{t) + 2ttW, 
where (3 = 1/T. The partition function is then written 
as a sum over the winding numbers, 



+ OC 



Z(q)= ^ lqW Zw, 



(4) 



W=-oo 



where Zw is a functional over the fields 4>w in the topo- 
logical class W given by 



D[4> w ]e 
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The action S[<f>] has the following form 



S[4>] = S d {4>] 



o 4£ c 



dr. 



(5) 



(6) 



The second term in the right hand side of Eq. (||) is 
the charging energy term. The first term describing 
the tunnel junction was obtained in Ref. ^ and can be 
conveniently written using the complex variable u = 
exp(27riTr) as follows, 

I dudu, Re (1 - 
Sd[<f>] = ~go <p 77T-T7 S 75 '-■ (7) 
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The integral in this equation is taken over the unit circle, 

\u\, \m\ = l. 



The minimum of the dissipative action (^) in the topo- 
logical sector W is equal to <7o|W|/2. In the trivial topo- 
logical sector, W = 0, it is achieved when 4>(t) — 0. In 
the sectors W = ±1 it is achieved on the instanton config- 
ur,ationsu which can be written in the complex notations 
asB 



exp(z0 2 ) = f(u) 



(8) 



where z (with \z\ < 1 for W = 1 and \z\ > 1 for W = — 1) 
is a complex number characterizing the position and the 
width of the instanton. In the topological sector with 
\W\ > 1 the dissipative action (0) is minimized on the 
multi-instanton configuration of the field <f> given by the 
product of \W\ single-instanton terms, as in the right 
hand side of Eq. (§). 

The topological sector with W = in Eq. does not 
contribute to the oscillatory part of the thermodynamic 
potential, f2(g) = —TlnZ(q). At relatively high tem- 
peratures, when the renormalized conductance is large, 
g> 1, the main non-zero contribution to the oscillatory 
part of comes from the terms with winding num- 
bers W = ±1 in Eq. (^). All other terms in Eq. (||) are 
exponentially small in g, and the renormalized charging 
energy in Eq. Q) can be expressed as 



E C (T) = 2T 



Zx 
Zo 



(9) 



We therefore concentrate below on the topological sector 
W = 1. 



III. SINGLE INSTANTON APPROXIMATION 

Since the value of the dissipative action on the instan- 
ton configuration <\> z in Eq. (gj) is independent of the in- 
stanton parameter z it is convenient to write the fields in 
the topological sector W = 1 in the form 



(10) 



where (j) z are massive fluctuations which are orthogonal 
to the two zero modes of the dissipative action, d6 z /dz 
and dcf> z /dz* . The renormalized charging energy (J9|) can 
be written as 



E C (T) = 2T 



d 2 z 



;Zi(z), 



(11) 



where Z\(z) is given by the following ratio of functional 
integrals over the massive modes only, 



Zi{z) 



J D[4> z ](l - \z\ 2 )J(z, 4> z ) exp(-5[^ + 4> z \) 



fDlfo} exp(-5[0 o ]) 



(12) 



where J(z,<f> z ) the Jacobian of the variable transforma- 
tion, Eq. @. 
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Below we compute th e ren ormalized charging energy 
Ec, Eq. (pi]). In section [II A we evaluate the functional 



integral over the massive mode s (|12|) in the Gaussian ap- 



proximation. In section 



[II B 



we obtain corrections to 



the Gaussian approximation in leadin g or der of pertur- 
bation theory in 1/go- Then, in section [TV] we employ the 
renormalization group to obtain Eq beyond the regime 
of validity of perturbation theory. 



A. Gaussian approximation 

To leading order in 1/go we may evaluate the ratio 
of the functional integrals in Eq. ( p"2| ) in the Gaussian 
approximation. To this end we expand the actions in the 
numerator and in the denominator in Eq. ( JT^ ) to second 
order in (f> and 4>o respectively. 

In the Matsubara basis, 



fluctuations of the massive modes, (p n are small as 1/go, 
see Eq. (|l6|). Therefore, for large g we can expand the 
determinant of the metric tensor in the powers of ip n . In 
leading order in 1/go, i.e. in the Gaussian approximation 
it is sufficient to evaluate detA(z,(/> z ) on the instanton 
configuration, (p n — 0. We show in Appendix [A| that 



J(z,0) = -\/detA(z,0) 



1 



f 



M 5 



(18) 



The quadratic form of the charging part of the action 
given by the second term in Eq. (p|) is not diagonal in the 
basis ([H]). However, for the purpose of evaluating the 
functional integral in the numerator of Eq. (12) one can 
neglect its off-diagonal elements. We show in Appendix 
[b| that their contribution is small as 1/go- The diagonal 
part of the charging action is given by Eq. (B6). Then, 
the ratio of the functional integrals in ( p^ ) reduces to the 
product 



E 



ip n U 



(13) 



the action in the denominator acquires the following di- 
agonal form, 



=3o^(n + an 2 )|^ n | 2 , 



(14) 



where we introduced the notation a — 27r F T . 

90 c 

To evaluate the functional integral in the numerator in 
Eq. (^) it is convenient to expand the massive fluctua- 
tions 4> z using the basis of Ref. || which in the complex 
notationscJ can be written as, 

4> z = Y / f'nU n f(u) + Y / ^nU n r(u). (15) 



n>0 



n<0 



Here f(u) is defined in Eq. (||). In this basis to second 
order in 4> z the dissipative part of the action, Eq. (Q) , has 
the following diagonal form, 



S': 



go 

2 



g ^2n\ip 



^ |2 



(16) 



The superscript i here refers to the presence of the in- 
stanton, whereas the superscript in Eq. ( |l4| ) denotes 
the trivial topological sector, W = 0. 

Instead of calculating the Jacobian J(z, (f> z ) in Eq. ( |l2] ) 
directly, we can express the integration measure through 
the metric tensor A(z, <j) z ) using the identitjilj 



detA(z, </> z ) JJd<£ r , 



(17) 



The metric tensor A(z,(/) z ) is presented in Appendix |A|. 
Fortunately we do not need to evaluate its determinant 
for an arbitrary configuration of the field <j) z . Indeed, the 



go 

7T 



n 

n=2 



an 



1 



2a\z\ 2 
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7T 



(19) 



where a was defined below Eq. ([14). Since the character- 
istic instanton frequencies are T/[l — \z\ 2 ) < Ec <C goEc 
the second term in the exponent of thi s ex pression can be 
neglected. Using the expression Eq. (B2) for the charg- 
ing action on the instanton configuration we obtain in 
the Gaussian approximation, 



2n 3 T 



ex P 77 



fjo 
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^r[i + M 2 ] 

Ec[l-\z\>] 



(20) 



Substituting this expression into Eq. (|Tl]) we obtain with 
logarithmic accuracy, 



E C (T) 



glEc 



In 



'T ' 






T 


exp | 


:-D 



where we introduced the notation 

E c 



To = 



2tt 2 



(21) 



(22) 



Equation (gTJ) coincides with the result of Ref. One 
can neglect the interaction between instantons arid ob- 
tain Eq. ( |2l| ) only when Z 1 (z) -C 1. In section IIIB we 
show that even in this regime the non-Gaussian correc- 
tions to the functional integral in Eq. ([l2]) lead to large 
corrections to the preexponential factor in Eq. (Ell) . 



B. Corrections to the Gaussian approximation 



The Gaussian approximation discussed in section [II A 
is asymptotically exact at go ~ * 00 ■ The corrections to it 
are small in 1/go and may be evaluated perturbatively. 
To obtain the leading 1/go correction it is sufficient to 
expand the Jacobian in Eq. ( |l2] ) to second order in 4> z 
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CMDO 



accuracy, 



FIG. 1: Diagrams representing 1/go corrections to j/ \ in 

Eq. (E5h . The vertices are proportional to go, and the Gaus- 
sian propagators of fields (j> depicted by the solid lines are 
proportional to 1/go- 



and to expand the actions in the numerator and in the 
denominator in Eq. ( fl2| ) to fourth order in <f> z and <j>o 
respectively: 



Sd[4n] = Si 2 + Si 3 + Si 



(23a) 
(23b) 



The terms S^ 2 and S^ 2 were defined in Eqs. ( fL4|) and 
( pf ) respectively, and 5 ~ 3 , S* ~ 4 , and are given by the 
following equations; 



cO _ 

— 



9o 
24 

9o 
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24 



Re 



Im 



dudui 
dudu\ 

WF 

dudu\ (<fi — 0i 
[2-ki) 2 (u — iti 



/(") 



o 4 



(it — Wi) 2 ' 

(0-^l) 3 

(w - Ul) 2 ' 



(24a) 
(24b) 
(24c) 



Here the function f(ui) was defined in Eq. (^), and we 
introduced the short hand notations (pi = 4>{ui). Substi- 
tuting Eqs. (H|) into Eq. (|l|) and using Eqs. (|l7|) and 
( IHf ) we obtain up to terms of order If go, 



Zi(z) 
Z?(z) 



(l-\z\')(^JdetA(z,4> 



where (. . .)i and (. . .)o denote averaging with respect to 
the Gaussian actions in the presence and in the absence 
of the instanton respectively. 

The first term in the right hand side of Eq. ( |25| ) needs 
to be evaluated to second order in <f> z . This is carried out 
in Appendix [A|. The result is given by Eq. (A4). The 
calculation of the other terms is facilitated by the use of 
Wick's theorem and reduces to evaluating the diagrams 
in Fig. [I]. The calculations are straightforward and are 
presented in Appendix Only fluctuations of <f> with 
relatively low frequencies, < T/(\ — |z| 2 ), contribute to 
the corresponding functional integrals. We can therefore 
neglect the charging energy term in the Gaussian action 
for the massive modes. The calculations are further sim- 
plified by observing that the two-point correlation func- 
tions in the diagrams in Fig. [l] become diagonal in the 
bases (|l3|) and ([L5|). We then obtain with logarithmic 



Z 1 {z) = Zf{z) 



go i go 



, (26) 



where T was defined in Eq. (p2|). Substituting Eqs. ( |2q ) 
and (H|) into ( |TT| ) and using Eq. (A4) we obtain for the 
renormalized charging energy, 



E C {T) 



9_ZEc ln n 

7T 2 T 



1 



ln^ 



In 2a 



-So/2 



5o 2gio 

(27) 

Equations (|26| ) and (p7|) represent the main results of 
this section. The first term in the right hand side of 
Eq. (j27]) coincides with the result of Ref . |[ and the others 
represent a perturbative correction arising from taking 
into account non-Gaussian fluctuations. Equation ( p7j ) is 
valid at relatively high temperatures, In Sr <C 50, when 
the non-Gaussian correction is small. We consider the 
region of lower temperatures, In ^ ~ g Q , in section ^ 
using the renormalization group approach. 



IV. 



RENORMALIZATION GROUP 



At low temperatures, when In becomes of the or- 
der of g , the non-Gaussian correction in Eq. ( |26| ) be- 
comes significant, and the perturbative expressions (|26| ) 
and ( pTj ) arc no longer valid. There is a wide temper- 
ature range in which the perturbative approach used in 
section |III fails but the renormalized conductance is still 

^> 1, and therefore the single 



5o 



21n^ 



instanton approximation, Eq. ( pd| ) is still valid. Below we 
apply the renormalization group to study the amplitude 
of the Coulomb blockade oscillations in this regime. 

From Eq. ( |20| ) we observe that Z\(z) depends ex- 
ponentially on the renormalized conductance, Z\{z) oc 
exp (— g(T)/2). Therefore to obtain the preexponential 
factor in the renormalized charging energy in Eq. jll| ) 
it is necessary to compute the renormalized conductance 
using the two-loop renormalization group. 

We obtain and solve t he tw o-loop RG equa tions for 
the conductance in Sec. IV A. Then, in Sec. 




obtain the renormalized charging energy Ec, Eq. 
by evaluating the functional integral Z\{z) in Eq 
with the aid of the renormalization group. 



A. Renormalized conductance 

To obtain the two-loop RG equations we expand the 
dissipative action in Eq. (^) to sixth order in (jjQ. Next 
we introduce a running frequency scale v and write 
4>o = 4>h + 4>i- Here 0q represents fast degrees of free- 
dom with Matsubara frequencies v n = Tn satisfying the 
inequality v <v n < 50Tb, and <Pq represents slow degrees 
of freedom with v n < v. The calculation of the renormal- 
ized conductance at frequency lu amounts to evaluating 
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FIG. 2: Diagrams of order 1/go for the renormalized conduc- 
tance g(v)- 



the diagrams in Fig. |], where internal lines correspond 
to the propagators of the fast degrees of freedom. 

The role of the charging term in the Gaussian action, 
Eq. (|l4|), amounts merely to the ultraviolet cutoff of the 
frequency integrals over the internal lines at frequencies 
goTo, where To was defined in Eq. (p2|). As a result we 
obtain, 

3( J/ )=ffo-21n In . (28) 

v go v 

From Eq. (|2^) we obtain the following two-loop renormal- 
ization group equation for the conductance in agreement 
with Ref. I 



dg{v) 



9(")' 



e=i n : 



A 



(29) 



where A = goT is an ultraviolet cutoff. The first term in 
the right hand side of Eq. (||^) reproduces the one-loop 
renormalization group equation of Refs. p|Ji"3|. 

By integrating the two-loop renormalization group 
equation, Eq. (29) from v = goTg to v = T we obtain 



Eq. (3b) which implicitly determines the temperature de- 
pendence of the renormalized conductance g(T). 



B. Renormalized charging energy 

Next, we apply the renormalization group to evalu- 
ate Zi(z) in Eq. (|l2|). We start by considering the 
case of the so-called line instantons, z = 0. Accord- 
ing to the RG procedure the functional integration is 
performed by separating 4> into fast and slow degrees 
of freedom, </> = $ + cf> s , where <f>f includes fluctua- 
tions with Matsubara frequencies v n in a narrow inter- 
val, v — Ae~£ < v n < A . Upon the integration over 
the fast degrees of freedom the remaining functional in- 
tegral acquires a factor resulting from integrating out 
the fast fluctuations of </>. This factor depends on the 
running coupling constant giy). The remaining action 
for the slow fluctuations is characterized by a renormal- 
ized conductance. Due to this multiplicative nature of 
the functional integration we can express the logarithmic 
derivative of Zi(0) through a function of the renormal- 
ized conductance g(v) only, dln ^°) — f[g(y)]. From 
the perturbative results (^6|) and ( |20| ) we obtain for the 
logarithmic derivative, 



dlnZi(O) 



1 



1 



9W 



(30) 



Using the two-loop renormalization group for the con- 
ductance, Eq. (P9[) we obtain, 



rflnZi(O) 

dgiT) 



1 



(31) 



The solution of this equation which satisfies the high tem- 
perature asymptotic (12^), (pp|) is given by, 



l(Q) = 2tt 3 P 



g(X) 



(32) 



One can directly check this expression against the pertur- 
bative result Eq. (|26|). To obtain Zi(0) to order 1/go we 
need to evaluate g(T) in the exponent of this equation to 
the two-loop order, and g(T) in the preexponential factor 
to the one-loop order only. 

To evaluate Z\(z) at finite z it is convenient to express 
it in the form 



Z t (z) = Z 1 (0)X{z)exp 



2 T[1 



Ec[l-\z\ 2 ] 



(33) 



where the charging energy on the instanton configuration 
is written out explicitly, cf. Eq. (|2(i|). 

The perturbative result ( p6| ) implies X{z) = 1 — 
^■ln(l — |z| 2 ), where the logarithmic term arises from 
the integration over the fluctuations of <fi with frequen- 
cies ranging between T and the instanton frequency 
v z = T/(l — N 2 ), see the discussions above Eq. ( p6| ) and 
in Appendix |Cj. Therefore in the renormalization group 
treatment the conductance go in the logarithmic correc- 
tion should be understood as the renormalized conduc- 
tance g{v z ) ~ g{T) + 21n-^=- at the instanton freque ncy 
v z . Using the considerations presented above Eq. ( |32] ) 
and the perturbative result, Eq. (^) we express the loga- 
rithmic derivative of X(z) through the renormalized con- 
ductance g{v z ) as, 



dhx[X(z)] 



1 



din 



g{v z ) 



(34) 



Using the one-loop renormalization group equation 
g7^7^ = 2 and the boundary condition X(0) = 1 we 
find, 



X(z) = 



g(T) 



(35) 



malized charging energy as 
d 2 z 



Substituting Eq. (33) into Eq. ( pT| ) we write the renor 

-X(z)exp 



E C {T) = 2TZ 1 (0) 



1-2 



* 2 T[l + \z\ 2 ] 



Ec[l-\z\ 2 ] 



Next, we express the integration measure through the 
differential of the renormalized conductance, dg{y z ), ac- 
cording to the one-loop RG equation, jAfp = — ^d ln(l — 
\z\ 2 ) = %dg(v z ). The resulting integral over the dg(ui z ) 
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ranges from qJuj z ) = gjT) to g(v z ) = g(T ) = g - 2 In go- 
Using Eqs. ( p2| ) and ( |3q ) we obtain for the renormalized 
charging energy, 



E C (T) 



Sir 2 



.,3/2 



(To 



7 3/2 



(T) 



-9(T)/2 



(36) 



where To was defined in Eq. (22). 

Next we express the temperature T through the renor- 
malized conductance g(T) using the solution of the two- 
loop RG equation, Eq. (BH). Substituting this result into 
Eq. ( fjif ) we obtain Eq. d3a|). 

Equations (^|) represent the main result of this pa- 
per. They determine the temperature dependence of 
the amplitude of Coulomb blockade oscillations in a wide 
temperature regime where the renormalized conductance 
g(T) is large, and the non-interacting instanton gas ap- 
proximations is valid. 

At relatively high temperatures, when In ^jfe <C go, 



Eq. (|3j) reproduces the perturbative result, Eq. (p7|). At 
lower temperatures, when In ^jP- ~ go, the renormalized 
charging energy in Eq. ( pa| ) significantly exceeds the pre- 
diction of the Gaussian approximation^, Eq. (|2^) . 

The non-interacting instanton gas approximation is 
valid while Zi/Z < 1. Let us extrapolate the results (||) 
to the low temperature limit, when Z\jZo ~ 1. This hap- 
pens at T ~ E c goe~ 9o/2 wheno(T) « 41n ff0 - Substitut- 
ing this conductance into Eq. (pah we obtain the following 
estimate for the amplitude of Coulomb blockade oscilla- 
tions in the ground state energy, Ec ~ •Ec l <7o e_So/ ' 2 - This 
exceeds the estimate arising from the Gaussian treatment 
of the fluctuations^! by a large factor, go- 



V. DISCUSSION 
A. Summary of the results 

We studied the amplitude, Ec(T), of the Coulomb 
blockade oscillations in the thermodynamic potential for 
a metallic grain connected to a lead by a tunneling con- 
tact with a large tunneling conductance-go- The effects of 
a finite mean level spacing in the grainEJ were neglected. 
We worked withing the non-interacting instanton gas ap- 
proximation. By applying perturbation theory we ob=. 
tained the leading 1/go correction to the Gaussian resultcl 
for the renormalized charging energy, Eq. ( p7| ) . Combin- 
ing the instanton analysis with the two-loop renormaliza- 
tion group we found the temperature dependence of the 
renormalized charging energy Ec{T), and of the renor- 
malized conductance g(T), Eq. (||), in a wide temperature 
range Ecgte~ ga/2 <T < E c - The use of the two-loop 
RG enables us to determine the preexponential factor in 
the temperature dependence of Ec (T) even at relatively 
low temperatures, when In ~ go- In this regime the 
renormalized charging energy in Eq. (|3a]) is significantly 
greater than the result of the Gaussian approximation, 

Eq. m>. 



Assuming that the amplitude of the Coulomb block- 
ade oscillations weakly depends on temperature at T < 
Ecgte' 9aj2 fromEq. (|) we obtain the following estimate 
for the amplitude of the oscillations in the ground state 
energy, Ec ~ -Ec<?o e ~ 90 ^ 2 ■ This estimateJs greater than 
the result of the Gaussian approximationO by a factor of 
go- In Ref. [?| the result for the integral over the massive 
fluctuations around the instanton configuration differed 
from that in Ref. || and in Eq. (|2C]). This lead to a dif- 
ferent estimate for the amplitude of the Coulomb block- 
ade oscillations in the ground state energy, Ecgo e90 ^ 2 ■ 
In Ref. |5| the instanton effects were ignored, and the 
estimate Ecgo^~ 90 ^ 2 f° r the amplitude of the ground 
state energy oscillations was obtained using the two-loop 
renormalization group. The topological effects, on the 
other hand, may drastically alter the results of a per- 
turbative RG consideration. This is well known, for ex- 
ample, ia the theory of antiferromagnctic Hciscnbcrg spin 
chains,E3 where the perturbative RG predicts the appear- 
ance of a spin gap for an arbitrary value of the spin S. 
However, for half-integer spin the topological effects are 
expected to lead to a gapless state.E3 



Applicability of the results and comparison with 
numerical studies 



Within the framework of the model considered here the 
thermodynamics of the grain is described by Eqs. (Q)- 
(0). The partition function, Eq. (^), depends on two 
dimensionless parameters, go and T/Ec- This. mode] has 
been studied numerically by several groups. EZIEj 

To compare our results with the numerical studies it 
is important to keep in mind the approximations that 
were made. These approximations determine the region 
of applicability of the instanton approach in the space of 
parameters go and T/ Ec- Below we review the approxi- 
mations made: 

i) The instanton configuration was found by neglecting 
the charging part of the action in comparison with the 

2tt 2 T 



dissipative part. This is justified if a = = <C I . 



T 
So T 

ii) To obtain the Gaussian result, Eq. (|2l|), we assumed 
that the charging action on the instanton configuration, 
second term in the exponent in Eq. ( p0| ) , is negligible for 
long instantons, z — > 0. This is valid for T <C Tq. In 
addition, Eq. (|l]) holds with logarithmic accuracy. This 
implies not only that ^ > 1 but that In ^ > 1. The 
low temperature regime is difficult to study numerically. 
At the lowest available temperature in Ref. |l7], In ^ = 
3.2. Thus, even at the lowest temperature, the result of 
Ref. f?j is not parametrically different from that of Ref. ||, 

Eq :.@- 

iii) The use of the single instanton approximation im- 
poses the lower bound, T* , on the temperature range of 
applicability of the results. In Ref. |h] the two-instanton 
contribution to the partition function was taken into ac- 
count in the one-loop approximation. It was shown that 
the single instanton approximation, Eq. (^|), applies as 
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long as < 1, i.e. for T* > E C (T*). It is easy to 
see from Eq. ( pl| ) that T* decreases as the dimension- 
less conductance go grows. The largest conductance for 
which the temperature dependent data are available in 

Rcf. is .go = 2tt 2 . Setting Ec ^T ] = 0.3 and using 
Eq. (|2l]) we find that JJ- w 10~ 2 . The temperature de- 
pendence of = 2n 2 Ec in Ref. [ll] saturates roughly at 
the same temperature. At this temperature lnSj- rs 1.6 
and quantitative comparison of our results with numeri- 
cal results is not justified. Qualitatively however, at this 
temperature the results of Eq. ( pl| ) and of Ref. are not 
very different. 

Thus quantitative comparison of numerical results 
with the results of different analytical approaches, 
Refs. 0J| and Eq. requires numerical studies at much 
lower temperatures and larger conductances than those 
available at present. 
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APPENDIX A: EVALUATION OF detA(z, 



To obtain the metric tensor A(z,<f> z ) in Eq. ( |17| ) we 
express the variation of the field <j> through the variables 
z and {(fin} with the aid of Eqs. ( |To| ) and (|l5|), 

ci d't'z x v"^ - n df s—^ . n df* 
dtp = ——dz+) (f n u —dz+y <p n u -^-oz 
oz oz ' oz 

n>0 n<0 



-T—dz + ytp n u n —Sz + ]>(p n u n —6z 
oz* ' oz* ' oz* 

n>0 n<0 
J2S<PnU n f(u) + J2 S VnU n f*(u). (Al) 



n>0 



n<0 



Below we omit the arguments z and 4> z of the metric 
tensor A. Its matrix elements are obtained from the fol- 
lowing relation, 

n,m 

+ A* zz Sz*Sz* + A zz 5z6z 

+ 2 J2 i A zm5z + A z * m 5z*} 6(p m .(A2) 
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A 
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\ 


A* 


A* 
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'•./ 



= A(z,5) 



FIG. 3: Matrix elements of the metric tensor A(z,<j> z ) in co- 
ordinates z, {<fin}- 



Substituting Eq. (Al) into Eq. (A2) it is straightforward 
to find all the elements of the metric tensor A(z,4> z ). 
Most of the matrix elements of A(z,<fi z ) vanish. The 
schematic form of this matrix is shown in Fig. |[ 

The determinant of A(z,(j) z ) in Eq. ([l7|) needs to be 
evaluated up to terms of order I /go- Therefore it is suf- 
ficient to find the matrix element A zz up to the li nea r 
order in (p. Integrating over the variable u in Eq. (A2) 
we obtain the following expressions for the elements of 
matrix A(z, <p z ), 



A z 



~n—l _ — n — 1 

^" l-lzl 2 ~ % '^- 

n>0 1 1 



\z\ 2 '- 

n<0 1 1 



(A3a) 



A, 



A — 



rr^p( 1+2 Ei^i 2 )' ( A3b ) 



E <P*n 
n>0 

- E ft 

n<Q 



^* v m— n— lz 



(m — n), to > 0, 
[m — n), m < 0, 



(A3c) 



Y *\n—m — l/ 



{n — m), m > 0, 



A z *m — 



- E K(z*Y 

n>0 

J2 {p* n {z*) n - m - x 6{n - m), to < 0. 

n<0 

(A3d) 

Here the ^-function should be understood as 6(0) = 0. 
We note that even for very long instanton, z — ► 0, some 
off-diagonal elements of matrix A(z, cf> z ) remain finite, see 

To evaluate the renormalized char ging energy Ec in 
the Gaussian approximation, see Sec. HI A, the Jacobian 
in Eq. (|l7]) may be evaluated on the instanton trajectory, 
4> z = 0. In this case the determinant of A(z, 4> z ) is readily 
evaluated, and we obtain Eq. (|lg|). 

In order to obtain the leading l/(?o correction to the 
Gaussian result for Ec, we need to calculate the square 
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root of the determinant of A(z, 4> z ) in Eq. ( B5h to sec- 
ond order in <j) z . To this end we expand det A(z,<p) to 
second order in its off-diagonal elements and perform 
the averaging with respect to the Gaussian fluctuations 
around the instanton, in Eq. (p5|). We show in 

Appendix [b] that the off-diagonal elements of the charg- 
ing part of the action may be neglected, and therefore 



(<Pn<P* r , 



= 5„r 



9o\n\ 



2tt 2 T 
Ec 



Using this result and Eq. (|A3|) we obtain 



(W + l) 2 



2\z\- 



detA(z, . 



1 



1 , 1 1 go r o 
. 1 + — In— — 

A \ go t 



(A4) 



APPENDIX B: 



CHARGING PART OF THE 
ACTION 



To obtain the charging energy, Eq. ( p7| ) we evalu- 
ated the functional integral over the massive fluctuations 
about the instanton in the basis of Eq. ( fl5| ) . We neglected 
the off-diagonal elements in the charging part of the ac- 
tion. In this appendix we show that this approximation 
is justified. Taking the off-diagonal matrix elements into 
account leads to a correction to Ec which is smaller than 
the non-Gaussian correction given by the last two terms 
in Eq. ( |27| ) by a factor of In ^S- . Therefore it can be 
neglected at low temperatures. 

To demonstrate this we rewrite the charging part of 
the action using the complex variable u = exp(27riTr) 
and Eq. @, 



^dr = ^l 



udu 




(Bl) 



was defined in 



Here the instanton configuration 

Eq. (||). First term in the right hand side of Eq. (|Bl| ) 



describes the charging part of the action on the instan- 
ton configuration and the last two terms describe small 
fluctuations around the instanton. Setting <j) z = and 



integrating over the variable u in Eq. (Bl) we obtain for 
the charging action on the instanton configuration, 



4E C E c \l-\z\ 2 



(B2) 



To calculate the second and the third terms in the right 
hand side of Eq. ( |Bl| ) we expand the massive fluctuations 
(f) z using Eq. (fL5|). Integrating over the variable u in 



Eq. (Bl) we obtain the following results 

niT f dd>,. <9cL . 2w 2 iT 



Er 



du du 



■udu — 



E c (l-\z\*) 



\n>0 n<0 J 



(B3) 



udu 



t: 2 T 



2E C J \du J E c 

^ [ftnfnm( z )f>-m+ t Pmfmn( z )<P-7 



(B4) 



n,m>0 



Here the function f nm (z) is given by the following ex- 
pression 

+(m + n) [z n -" l d{n - m + 1) + {z*) m ~ n 6{m - n)] 
l + |z| 2 " 



+z 



n — m + 



l-UI 



+(z*y 



m — n + 



1 



1- z 



(n-m + 2) 



(m-n-1). (B5) 



Here again 9(0) = 0. From Eqs. (|B4|) and flBJ) it fol- 
lows that the diagonal part of the quadratic form of the 
charging action is given by the following expression, 



S- 



diag 



2tt 2 T 
~Ec7 



S |2 



ri>0 



(n + l) 2 + 



2\z[ 



1- z 



(B6) 



The result of Eq. @ was used in Eq. (jT|). 

We now consider the contribution to the renormalized 
charging energy Ec(T) in Eq. ( p7| ) from the term S2 in 
Eq. (B3). We expand the numerator in the right hand 
side of Eq. ( |l2| ) to second order in £2. As a result we 
obtain the following correction to the right hand side of 
Eq. ©, 



<^ G (*) = #f^ex P 



, g0 ^ 2 T[1 + |z| 2 ] 
2 E c [l-\z\ 2 ] 



2ir 3 T 2 

J "(B7) 

Here the angular brackets, (. . denote averaging with 
respect to the Gaussian action, Eq. Jl6|). From Eq. (^ ) 
we obtain, 



2^ 2 T 



E c [l-\z\< 



Tl>0 



Using the fact that (l^ril 2 )^ = (<7o 1^1) 1 and substitut- 
ing Eqs. (B7), (B8) into Eq. ( |IT1) after integration over 
the z we obtain the following correction SEc(T) to the 
renormalized charging energy Ec(T) in Eq. (B7 



5E C (T) = ^ 



lnf- 
9o 



exp 



(-?)■ < B9 > 
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This correction is proportional to the first power of the 
large logarithm ln(T /T). Therefore it is much smaller 
than the second term in the right hand side of Eq. ( pjj ) 
and can be neglected. It is straightforward to check 
that the contribution to the renormalizcd charging en- 
ergy Ec(T) in Eq. @ from the term S 3 in Eq. (gJHs 
small as 1/qq in comparison with the result of Eq. (B9) 
and therefore can also be neglected. 



As is verified by the subsequent calculations, only fluc- 
tuations of <j) with frequencies below T/(\ — \z\ 2 ) con- 
tribute to the corresponding functional integrals. There- 
fore we neglect the charging energy term in the Gaussian 
action for the massive fluctuations of cf>. It is then conve- 
nient to perform the calculations in the bases ([l3|), (|l5| ) 
which diagonalizeO the dissipative action. For example, 
for the correlator ((</) — </>i)(</> 2 — 4>3))i we obtain the fol- 
lowing expression, 



APPENDIX C: EVALUATION OF THE LAST 



THREE TERMS IN EQ. (25) 



We employ Wick's theorem to rewrite the last three 
averages in Eq. (E9) in the following form, 



l^nl 2 )* 



n>l 



xRe 



— — — -h n (u,ui)h n (u 2 ,U3) 



(CI) 



ii) 3 (02-03) a ). = 6((0-0i)(0 2 -0 3 ; 



+9((0~0 1 )(0 2 -03)).((0-0l) 2 ).((^2-03) 2 )., 



((0-0 1 ) 4 ) o = 3((0-^ 1 ) 2 ) 2 . 

Different pairings in this equation can be represented by 
diagrams in Fig. El 



Here we introduced the notation 



h n (u,ui) = 1 



u 

Ui 



m 



(C2) 



All other correlation functions may be expressed through 
(\{pn\ 2 )i and (|ty? n | 2 )o hi a similar way. Since upon neglect- 
ing the charging action {\<f> n \ 2 )i = (|<^n| 2 )o = (50M) -1 
we can write the averages in the right hand side of 
Eq. (Eq) in the following form: 



at 



< s ».) 2 >, 



1 y- _L 

i,n\ >1 
00 1 

E - 

Z ^ 7777 



dudu\ 



dudu\ 



-Re 



On ' — ' nni / (2vTi) 2 (u — ui ) : 

MU n,m> 2 1 17 v ' v i; 



-Re 



m 

/(«l) 



Re [/i„(it, ui)] Re [/i ni (u, ui)] , 



It 


n- 






u 


m- 






Re 


1 - 














Ui 





(C3a) 
(C3b) 



2 1 / duduidu^duz 

3g a " nnina / (27r) 4 (u - ui) 2 (u 2 - U3) 2 

n,ni,n2>l / 



Im 



/(«) 



: (6Re[/i ni (u,ui)]Re[/i n Ju 2 ,U3)]Re[y„] + Re[y„]Re[y ni ]Re 



(C3c) 



where y nj denotes the following function, 
"* /(ui) 



Dm 



Ml 
"2 



/(«a) 



h ni (u,ui)h ni (u 2 ,u 3 ). 



(C4) 



The two terms in the last line of Eq. ( C3c ) correspond 
to the second and third diagrams in Fig. [1] respectively. 

Next we perform the integration over Ui in Eq. (|C3| ). 
The first term in the last line of Eq. (C3c) vanishes. The 
remaining integrals lead to the results: 



— y 

250 ^ 

n,ni >1 



2min(n, m) — |z| z 



(C5a) 



212 00 



E (C5b) 



2 6# wrn«2 

n,ni ,ri2> 1 

x [2(rc 2 - n - n 1 ) 2 |«| 2(na ~ Bl-n - 1) e(n a - n - m) 
+(n - n 2 - n 1 ) 2 |z| 2 ("-" 1 -" 2 - 1 ^(n - n 2 - m)], 
1 min(n, ni) 



Wo 



E 



5o ram 

n,n± >1 



(C5c) 



Evaluating the sums in Eq. (C5) with logarithmic accu- 
racy and using Eq. (E3) we obtain Eq. (Eq). 
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